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Abstract 

In this work, we study Hawking radiation from a general static black hole due to tunnelling 
of particle having nonzero mass. Hawking temperature has been calculated using both the 
tunnelling method and the Hamilton- Jacobi method and the results do no depend on the 
mass of the particle. Due to complicated form of the equations involved quantum corrections 
are evaluated only upto first order and it is possible to intcrprctc the correction term as loop 
corrections of the back reaction effects in the space time. Finally, modified expressions 
for entropy and surface gravity has been evaluated and it is found that the leading order 
correction term to the entropy is not necessarily logarithmic in nature. 
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1 Introduction 

Classically, nothing can escape from a black hole across its event horizon. But in 1974, there was a drastic 
change in view when Hawking [T| [2] showed that black holes are not totally black, there is a continuous 
emission of radiation from a black hole identical to black body radiation [3J|3] having temperature T = 
K being the surface gravity of the black hole. Subsequently, there has been lot of attraction to this issue 
and various approaches has been developed to derive Hawking radiation and corresponding temperature 
[3J HI O |H1 [7] . In recent past a semiclassical approach [5J H] has been formulated considering Hawking 
radiation as tunnelling phenomena across the horizon. The basic idea behind the semiclassical approaches 
is to evaluate the imaginary part of the action (S) for the (classically forbidden) process of s-wave emission 
across the horizon and is then related to the Boltzmann's factor to obtain Hawking radiation. In fact, 
The studies of black hole radiation can be classified into two groups. The first approach developed by 
Parikh and Wilczek [5] is based on the heuristic pictures of visualization of the source of radiation as 
tunnelling and is known as radial null geodesic method. The essence of this method is to calculate the 
imaginary part of the action for the s-wave emission using the radial null geodesic equation. This method 
is limited to massless particles only and is applicable to such coordiate system only in which there is 
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no singularity across the horizon. The alternative way of looking into this aspect is known as complex 
paths method. Srinivasan et.al. [TQl HI] did the pioneering work for these Hamilton Jacobi(HJ) approach 
considering emitted particle without self-gravitation and assuming the action satisfies the HJ equation. 
Tunnelling of both massless and non-zero mass particles are possible in this approach and it is applicable 
to any coordinate system to describe the black hole. 

Most of the studies [121 E3] dealing with Hawking radiation are concentrated to semiclassical analysis 
and are confined to massless particles. Recently, Banerjee et. al. [2] initiated the calculation of Hawking 
temperature with quantum correction and have evaluated the modified Hawking temperature, surface 
gravity and entropy. Also they have interpreted the parameters involved in quantum corrections by loop 
corrections to backreaction effects in the space time. 

In the present work, we consider tunnelling of a nonzero mass particle using both the approaches 
in the semiclassical approach of Parikh et. al. 8 (which if so far applicable for massless particles only) 
we have evaluated the radial velocity of the wave front of the de Broglie's wave as spherical wave and 
evaluated the imaginary part of the action. The semiclassical Hawking temperature is then evaluated in 
section 2 by comparing with Boltzmann factor. Quantum correction to the Hawking temperature has 
been evaluated using HJ approach in section 3 and the semiclassical part of the temperature agrees with 
that in section 2. The evolution of the horizon is studied in section 4. We have shown that the static black 
hole in the tunnelling phase behaves like Vaidya type black hole and tunnelling of massless and massive 
particles have identical radial velocity. Finally, in section 5 we have given modified expressions for surface 
gravity and entropy with quantum corrections. Also the parameters are analysed from the point of view 
of loop corrections to back reaction effect of the space time. The paper ends with a summary in section 
6. 



2 Semiclassical Tunnelling approach 

The section deals with tunnelling of nonzero mass particle considering the picture of Hawking radiation. 
In fact the method correlates the imaginary part of the action for the classically forbidden process of de 
Broglie wave emission across the horizon with the Boltzmann factor for the black body radiation at the 
Hawking temperature. 

A general static black hole metric has the form 

ds 2 = -A(r)dt 2 + -^--dr 2 + r 2 d^l\ (1) 
v ' B(r) z v ; 

where the location of the horizon r — r/j is given by A(rh) = — B(rh). To remove the coordinate 
singularity at the horizon we consider the above metric in Painleve coordinates |15j as 



ds 2 = -A(r)dt 2 + 2 J -^y (1 - B{r))dtdr + dr 2 + r 2 dVi\ (2) 



which is obtained from equation ([1} by the transformation 



dt^dt-J ], B ^\ dr (3) 
V A(r)B(r) y ' 



The Painleve coordinates has a number of interesting features, namely 
(i) the spatial geometry of any i=constant hypersurface is flat, 
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(ii) the boundary geometry for any fixed radius is identical to that of the metric in equation ([T]) and 

(iii) the metric has time reversal asymmetry. 

Unlike the massless particle moving along the null geodesic a non-zero mass particle is treated as 
a de Broglie wave (s-wave). Further, to preserve the spherical symmetry of the space-time during the 
tunnelling process, we treat it as a spherical wave with velocity of the wave front as |16j 



• = = I Mr)B{r) 

2g 01 2\jl-B(r) 

Then acording to Parikh and Wilczek [8] the imaginary part of the action is obtained as 



(4) 



ImS — Im / p r dr — Im / / dp' r dr = Im / < / —, — > dr (5) 





' r H dH') 


n 






Jo % J 



where (r, p r ) are canonical pair and Hamilton's equation f — 



dp, 



has been used in the last step. 



In ordinary quantum mechanics, the action of a tunnelled particle in a potential barrier having energy 
larger than the energy of the particle will be imaginary as p r — ^J2m (E — V). For the present static 
black hole, the tunnelling of a nonzero particle of energy 1 uj' implies dH' integration over all possible 
values of energy of outgoing particle from zero to lj. Further, as we are dealing with tunnelling across 
the black hole horizon so using Taylor series expansion about the horizon Th we write 



Mr) - ^\ r=rh (r-r h )+0(r-r h f 



(G) 



So from (UJ) the velocity of the wave front in the neighbourhood of the horizon can be approximated as 



dr _ y/A'(r h )B'(r h )(r-r h ) 

dt ~ 1 ~ [ ' 



lyJl-B'(r h )(r-r h ) + 0(r-r h y 



where dash denotes differentiation w.r.t. 'r'. 

Using (J7]) in ([S]) and evaluating the contour integration over V we obtain 

ImS= , 2 ™ (8) 
yJA'(r h )B'(r h ) 

where the contour for r-integration is chosen on the upper half complex plane to avoid the coordinate 
singularity at 77,. Thus the tunnelling probability is given by 

r~exp(-!/m4 =exp{ ^ ) (9) 

I » J \ hy/A'(r h )B'(r h ) J 

which in turn matching with the Boltzmann factor exp { — y}, we obtain the Hawking temperature as 



hy/A'(r h )B'(r h ) 

4^ (10) 

This is the welknown semiclassical Hawking temperature for tunnelling of a particle across the static 
black given by equation (JXJ) . Note that T# is independent of the mass of the particle. 
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3 Hamilton-Jacobi method : WKB approximation 



The HJ method will be employed to consider tunnelling of particle having nonzero mass beyond semi- 
classical approximation. We shall restrict ourselves to only loading order quantum correction due to 
complicated form of the equations involved. 



The Klein Gordan (KG) equation for a scalar field ip describing a scalar particle of mass mo has the 
form [10] 

,2 ~ 



(ii) 

where the d (box) operator is evaluated in the background of the black hole metric (Q}. Due to spherical 
symmetry, without any loss of generality we may restrict ourselves only to (t, r)- sector on the space-time, 
i.e., we fix our attention to two dimentional black hole problem. So the explicit form of the above KG 
equation is 



1 d 2 i) d 2 ip 



1 d 



dtp 



A dt 2 



dr 



2A dr 



Or 







(12) 



In W.K.B. approximation, we substitute the standard ansatz for the semiclassical wave-function namely, 



fl>(t, r) 



exp 



- h S(r, t) 



then the action S will satisfy the following differential equation 



1 (21 
A \~dt 



B 



dS 



dr 



1 d 2 S 



A dt 2 



B 



2 S 



dr 2 



1 d OS 
2ATr {AB) ^ 



= 



Now, quantum corrections are introduced by expanding S in powers of Planck constant 'Ti'as 

S(r, t) = S (r, i) + ^^ fe 5 fe (r, t) 



(13) 



(14) 



(15) 



where So is the semiclassical action and k is a positive integer. Using this ansatz for S into the differential 
equation (1141) and equating diferent powers of h in both sides we obtain the following set of partial 
differential equations (p.d.e.) : 



n 1 

I (21l 
a\ at 



2dS dS 1 o dS dSi. 

IB— h l 



A dt dt 



dr dr 

2 



- B 

1 d 2 S 
A dt 2 



dSo 
dr 



- = 



B 



d 2 Sn 



A dt dt \ dr 



or or 



dr 2 

1 d 2 s l 

A dt 2 



1 d (AR\ dS ° 
2Adr {AB) ~dr~ 



B 



d % S x 
dr 2 



2Adr {AB) ^r~ 



as. 



(16) 
(17) 

= 
(18) 



and so on. 

The differential equation (jTS]) for the semiclassical action Sq can be solved by choosing the separable 
ansatz 

S (r, t) =Lo Q t + D Q (r) (19) 
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where 



A>M = ± 



AB ° dr = ±I (say) 



(20) 



Here +/— sign corresponds to ingoing/outgoing scalar particle substituting this choice for So into equation 
(fT7|) . We have the differential equation for the first order correction Si as 



dSi I 2 ~ 8S ,lm§ 



/ cL4 
Wq — dr 



using similar ansatz for Si namely, 



Si{r, t)=u 1 t + D 1 (r) 



we have from equation (|21j) 

£>i(r) = ±uj uji J 
= ±h- 



dr 



* In 

y/AB^-m'A) 



■ In or — m£A\ 



(say) 



Due to complicated form, retaining terms up to first order quantum corrections, i.e. 

S = S + KSi = (u + hut) t + {D (r) + hD^r)} 



(21) 
(22) 

(23) 
(24) 



Thus the wave functions denoting ingoing and outgoing solutions of the KG equation (fT2"j) using (fT3f and 
(l24l are of the form 



exp 



ip ou t = exp 



i 

h 

i 

~h 



(ujo + Tiui ) t + Jo + Uh — — In 



Icj 2 — ml\A\ 



ih 

(cjq + Huji) t — Iq — hli — In | c<j 2 — TOqA| 



(25) 



(26) 



As during tunnelling across the horizon (r, t) changes their coordinate nature so to overcome this difficulty 
we take into consideration that the time coordinate has an imaginary part in crossing the horizon and 
accordingly the temporal part has contribution to the probabilities [141 117) . Thus the incoming and 
outgoing probabilities are given by 



Pit 



and 



\lpi 



Wo 



exp 



exp 



— ^ Im [(u>o + huii) t] + Im [la + hli] — Im 



ih ,|2 2 a\ 

— In (J — m n A\ 



Im [(loq + fiuji) t] — Im [Iq + hl\\ — Im 



ih 



hi 



\uj 2 ~ m^Al 



(27) 



(28) 



As in the classical limit the horizon behaves as a one way membrane, no barrier for incoming particles 
so in the limit, i.e., 

(29) 



lim Pi n = 1 



Hence from (l27l) we have 



Im (u>ot) = Im Iq 
and 

Im (uit + h) = Im [| In \lo 2 - m^A\] 



(30) 
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As a result P ou t turns out to be in simple form as 



Pout = exp {-jjm (J + Mi)} 
I I— 



exp 



exp 



4,ui 



-Im 



4toQ I 1 

n 1 1 



\1 



f"' 
Jo 



AS 



dr 



V 



AS 1- 



So from the "detailed balance" principle [TU1 EJ HH] namely 

, E~ 

P out = exp 

the temperature of the black hole is given by 



j I Pin = exp<j - — 



T ft £ 

Jb = -; 

47T Wq 



1+^1 



J y/A' (r h )B' (r h ) 



= 1 



H 



(31) 



(32) 



(33) 



(34) 



where energy of the radiated particle is taken to be u; (i- e -> E = ojq). Thus the temperature of radiation 
does not depend on the radiated particle at least upto the first order quantum correction and it agrees 
with the claim in the literature. 



4 Evolution of the horizon during tunnelling 

At the beginning we try to address the question : "Do static black holes evolve to Vaidya black hole 
during Hawking radiation?" 

To have a possible answer to this question we start with a general static spherically symmetric black 
hole metric 

ds 2 = -f{r)dt 2 + J^+ r 2 dn 2 2 (35) 



where f(r) without any loss of generality can be chosen as 

r 

so that at the horizon r — rh, f(fh) = an d we have 

2M(r h )=r h 



(36) 



(37) 



As before in Painleve coordinates (i.e., dt — >■ dt — ^—j^y^-dr) the metric (|3"5|) can be written in singularity 



free form as 



ds 2 = ~f{r)dt 2 + 2 >/! - f(r)dtdr + dr 2 + r 2 dfll 



(38) 
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Now to match with the Vaidya black hole we make another transformation similar to the ingoing 
Eddington-Finkelstein type coordinates (v, r, 0, </>) with 



/dr 
; (39) 

Accordingly, the metric ([55)1 becomes 

ds 2 = ~f(r)dv 2 + 2dvdr + dr 2 df% (40) 



If M in equation (|36|) is now a function of v then equation (|40l) . which is the line element of the static 
black hole, can be identified as the Vaidya type black hole metric. 

Basically, the tunnelling process described by Parikh and Wilczek [8] is the image of that originally 
described by Hawking. According to Hawking [I], a virtual particle pair is created just outside the 
horizon due vacuum fluctuations. The negative energy virtual particle, which is forbidden in the exterior 
to the black hole, tunnels to the interior while the positive energy particle materializes as a real particle 
and escape to infinity. On the other hand, Parikh and Wilczek considered pair creation inside the 
event horizon and the positive energy virtual particle tunnels out to the exterior. However in both the 
processes the positive energy virtual particle becomes real one outside the black hole and seems to be 
radiated (Hawking radiation) from the black hole while the negative energy particle is absorbed by the 
black hole. As a result, there is a decrease in the mass of the black hole and consequently, the black 
hole will shrink due to the effect of self-gravitation. The barrier in the tunnelling process is set due to 
contraction of the event horizon or the self gravity effect of the tunnelling particle and the two turning 
points can be identified as 

U = r hi and tj = r hf (41) 

where and r/j, are the positions of the event horizon before and after the tunnelling. 
For example in case of Schwarzschild black hole 

U = r hi = 2M and t f = r hj = 2(M - ui) (42) 

where u) is the mass (or energy) of the tunnel particle. 

Thus in course of tunnelling, the static black hole becomes a dynamic one and then they can be 
identified as Vaidya type black hole gven in equation (|4T))) . 

We shall now show that the radial velocity of the escaping particle does not depend on its mass- it is 
same for both massless or non-zero mass particle. A massless particle moves along a radial null geodesic 
which is characterized by ds 2 = dO = dip = 0, the velocity of the wave front is given by 

dr 1 

i = s'M < 43 > 

As a nonzero mass particle moves along a time-like geodesic so we can treat it as a de Broglie wave. 
Further, to maintain the spherical symmetry of the space time in course of the tunnelling process, we 
assume it as a spherical wave and consequently, the velocity of the wave front of the de Broglie s-wave is 

m 

dv 2 501 2 

Hence the velocities of the wave fronts for the massless or massive particles are identical. 
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Further, due to self gravitation, the mass M in f(r) given by equation (|36l) corresponding to the 
Vadya metric (|40j) will become M — uj where a particle with a mass of w tunnels out the horizon. Hence 
the rate of decrease of the horizon is given by 



Thus in the originsl position of the event horizon (i.e., r = r^) the shrinking velocity of the event horizon 



is 



fh= -^ m = -^ WF) =-- (46) 

But for Vaidya black hole given by equation (|40|) . the local event horizon (a null hypersurface, preserving 
the symmetry of the space time) is given by |19j 

r h = - ^ 47 

So using (|4T)]) the position of the event horizon (after tunnelling) is given by (neglecting square and higher 
powers of lo in comparison to M) 

r h = ^^~2(M-o;)=r(l-/(r))-2 W (48) 

M 

Thus we have t% — Th and tf = rh — 2lu are the two turning points of the tunnelling process and the 
contraction velocity of the event horizon is given by (146[) . 

5 Modified Surface gravity and entropy of the black hole 

In this section we derive expressions for surface gravity and entropy of the black hole considering to 
leading order quantum correction for the Hawking temperature. Using the standard relation between the 
surface gravity (k c ) and the Hawking temperature Th, i.e., 



Z7T 



the quantum corrected surface gravity is given by 



k q = 2ttT b = k c [ 1 + — ) (49) 



(50) 



Assuming 

uj\ e 
ujq M 

where e is a small dimensionless parameter then if we retain terms upto first order in e then from (|49[) 

K Q = «cU- j^J (51) 

According to ref. [131 HO] the above expression can be considered as the one loop back reaction effects in 
the space time. Also the dimensionless parameter e can be related to the trace anomaly 14 . 
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Further, if in the expansion in equation (|49[) we retain terms upto second order in e then we get 



^ = ^-W + M*) (52) 

This is nothing but the correlated form of the surface gravity of a blck hole due to two loop back reaction 
effects in the space time and e 2 is related to the contribution from the second loop. Therefore, expanding 
in higher order terms one can reproduce other higher order loop back reaction effects in the space time. 
The usual Bekenstein-Hawking area law 

S C = ± (53) 

where A, the area of the horizon, is valid at the semiclassical level. To find the modified form of the entropy 
when quantum corrections are taken into consideration, we take care of the first law of thermodynamics. 
Essentially, it is a energy conservation relation in which a change in the black hole mass M is related to 
the change in its entropy by the relation(assuming uncharged and nonrotating black hole) 

dM = T B dS B (54) 



i.e., 



dM 



dM 



\ u 



where Th is given by equation (If 01) . Writing explicitly we have 

Note that the second term on the r.h.s. of equation (|56p which is the leading order correction term can 
not always give rise to the logarithmic correction — it depends on the nature of the derivative of the 
metric coefficients at the horizon. Therefore, in general for a static black hole the leading order quantum 
correction may not give rise to a logarithmic form. 



6 Summary 

In this work, we consider tunnelling of a non zero mass particle across the horizon of a general static black 
hole. First of all, we have used the semiclassical tunnelling approach of Parikh and Wilczek to determine 
the Hawking temperature. Normally, this approach is used for zero-mass particle travelling along radial 
null geodesic and hence this method is also known as radial null geodesic approach. In section 2 we have 
used the tunnelling method for no zero mass particle where we have treated the particle as a de Broglie 
wave and the radial velocity of the wave front is evaluated considering it as a spherical wave. Then HJ 
method is used to incorporate quantum corrections. Due to complicated form of the equations involved 
we restrict ourselves to the leading order quantum corrections only. The results do not depend on the 
mass of the tunnelling particle and the semiclassical Hawking temperature agrees with the result of the 
previous section. Subsequently, we have shown that all static black holes with g o — 9rr > otherwise gen- 
eral become Vaidya-type black hole during tunnelling process and we have shown that the radial velocity 
of the wave front corresponding to particle of zero mass or non zero mass are identical. Finally, we have 
evaluated modified surface gravity and entropy of the black hole with quantum correction. We have seen 
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that for general static black hole first order quantum correction may not be logarithmic in nature. Also 
with the proper choice of the arbitrary parameters involved it is possible to show different order loop 
corrections of the back reaction effects in the space-time and the dimensionless parameter may be related 
to the trace anomaly. 
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